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CCD BV and 2MASS photometric study of the open 
cluster NGC 1513 
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Abstract We present CCD BV and JHK S 2MASS 
photometric data for the open cluster NGC 1513. We 
observed 609 stars in the direction of the cluster up 
to a limiting magnitude of V ~ 19 mag. The star 
count method shows that the centre of the cluster lies 
at a 20 oo = 04' l 09 m 36 s , 6 2 ooo = 49°28'43" and its angu- 
lar size is r = 10 arcmin. The optical and near-infrared 
two-colour diagrams reveal the colour excesses in the 
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direction of the cluster as E(B - V) = 0.68 ± 0.06, 
E(J - H) = 0.21 ± 0.02 and E(J - K s ) = 0.33 ± 0.04 
mag. These results are consistent with normal inter- 
stellar extinction values. Optical and near-infrared 
Zero Age Main-Sequences (ZAMS) provided an aver- 
age distance modulus of (m — M) = 10.80 ± 0.13 mag, 
which can be translated into a distance of 1440 ± 80 
pc. Finally, using Padova isochroncs wc determined the 
metallicity and age of the cluster as Z = 0.015 ± 0.004 
QM/H] = -0.10±0.10dex) and log{t/yr) = 8.40±0.04, 
respectively. 

Keywords Galaxy: Open Cluster and associations: 
individual: NGC 1513 stars: interstellar extinction 



1 Introduction 

Systematic studies of open clusters help understand the 
galactic structure and star formation processes as well 
as stellar evolution. By utilizing colour-magnitude dia- 
grams of the stars observed in the optical/near- infrared 
(NIR) bands, it is possible to determine the underlying 
properties of open clusters such as age, metal abun- 
dance and distance. In this study, we continue this ef- 
fort by analyzing the optical and NIR data of the open 
cluster NGC 1513. 

The young open cluster NGC 1513 = C 0406+493 
(a 2000 = 04' i 09 m 10 s , S 200 o = 49°3l'oo", I = 152°.50, 
b = — 1°.66) is classified as a Triimpler class II lm, 
a moderately rich cluster with a low central concen- 
tratio n. NGC 1 513 wa s first studied astrometrically 



by iBronnikoval (|1958a| ) who determined the proper 
motions using a single p air of plates with an epoch 



differe nce of 55 years. Barhatova fc Drjakhlushinal 
(|l960h published photographic and photovisual mag- 
nitudes of 49 stars from Bronnikova's (1958b) list. 
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Del Rio fc Huestamendi a (1988) obtained the first pho- 
toelectric UBV magnitudes of 31 stars and photo- 
graphic RGU magnitudes for 116 stars in the clus- 
ter region. They determined colour excess, dis- 
tance and the age of the cluster from RGU data as 
E(G - R) = 0.93 mag, 1320 pc and log(t/yr) = 
8.18, respectively. To test stellar evolution models, 
IFrandsen fc ArentoftJ dl998h studied NGC 1513 using 
BV photometry. Frolov et al. ( 2002J ) studied the astro- 
metric data-sets of 333 stars in the direction of the clus- 
ter, and showed that 33 of those stars are most proba- 
bly cluster members. Using BV photometric data, they 
also determined the cluster's metal abundance to be 
about solar metallicity and its ag e as \og(t/yr) — 8.40. 



Macieiewski fc Niedzielskil (120071 ) studied 42 open clus 



ters using CCD BV photometry and obtained the struc- 
tural and astrophysical parameters for NGC 1513. Ac- 
cording to their measurements the cluster's colour ex- 
cess, distance modulus and age are E(B — V) = 0.76 
mag, (m — M) = 12.96 mag and \og(t/yr) = 7.4, re- 
spectively. 

In this study we observed the open cluster NGC 1513 
with BV filters, and matched our results with 2MASS 
photometry. We determined the colour excesses in op- 
tical and near-infrared region. Then, we used optical 
and near-infrared data to obtain the cluster's distance 
modulus, metal abundance and age. 

The paper is organized as follows. We present the 
observations and data reductions in Section 2, while 
in Section 3, we describe the data analysis. Finally, 
Section 4 contains the conclusions of our study. 



2 Observations and Data Reduction 

2.1 Optical Data 

CCD BV photometric observations of NGC 1513 were 
made in 8th and 9th of October 2004 at the TUBITAK 
National Observatory (TUG) using the 1.5-m Russian- 
Turkish Telescope RTT150 and ANDOR DW436 CCD 
camera (back illuminated, 2kx2k pixels, 13.5x13.5 
/zm). The resulting image on the CCD has a field of 
view of 8' x 8'. Since the spatial diameter of NGC 1513 
is approximately 15' x 15' we divided the field into four 
equal subfields and created a mosaic image from these. 
Coordinates of the subfields and observation summary 
are given in Table 1 and the subfields are shown in Fig. 
1. 

For each subficld we obtained six images in B and V 
bands, which are then combined for each filter and me- 
dian images are used for further analysis. Obtained im- 
ages were reduced using the computing facilities avail- 



able at TUG, Turkey The standard IRAF0 routines 
were utilized for prereductio n, and the IRAF ve rsion 
of the DAOPHOT package (|Stetsonl 1 19871 . Il992h was 
used with a quadratically varying point spread function 
(PSF) to derive positions and magnitudes for the stars. 
To determine the PSF, we used several well-isolated 
stars in the entire frame. Output catalogues for each 
frame were aligned in position and magnitude, and final 
(instrumental) magnitudes were computed as weighted 
averages of the individual values. Magnitudes of the 
stars brighter than V = 10 could not be measured due 
to saturation of the detector pixels. 

The instrumental b and v magnitudes were trans- 
formed into standard Johnson B and V magnitudes us- 
ing fitting coefficients derived from observations of the 
standard field stars whose mag nitu des were given by 
)el Rio fc Huestamendial (|l988l ) and Frandsen fc ArentoftJ 
1998|) for photoelectric UBV and CCD UBV photome- 
try, and taking airmass correc tions into account. 15 and 
16 sta ndar d stars given by Del Rio fc Huestamendial 
( 19881 ) and IFrandsen fc ArentoftJ (|l998l) have a mag- 
nitudc and colour range of 12 < V < 18 and 0.55 < 
B — V < 1.70 mag. Errors of magnitude and colours 
for those stars are gi ven as ±0.02 and ±0.01 for V 
and B-V, respectively (IDel Rio fc Huestamendialll988 : 
Frandsen fc Arentoft 19981 ). Due to our subfield se- 
lection, we had 8-11 standard stars in each field. We 
then used these standard stars to estimate the standard 
magnitude of the other stars in the field. 

The m agnitude and colour differences betw een our 

values and Del Rio fc Huestamendia ( 1988J ) and Frandsen fc ArentoftJ 
( 1998J ) are shown in Fig. 2. The mean magnitude and 
colour differences and their respective standard devia- 
tions are 0.005, -0.011 and 0.07, 0.07 mag. 

The internal errors, as derived from DAOPHOT, in 
magnitude and colour are plotted against V magnitude 
in Fig. 3 and the mean values of the errors are listed 
in Table 2. Figure 3 shows that photometric error is 
< 0.015 mag at V ~ 19 mag while the colour error is 
0.024 mag. 

We compared our results with a rec ent BV photo- 

metric study on NGC 1513, performed by Macieiewski fc Niedzielski^ 
( 20071 ). 452 out of 609 stars in our sample matched 
with Maciejewski & Niedzielski's (2007). The mag- 
nitude and colour differences between our results and 
those of Maciejewski & Niedzielskis (2007) are given in 
Fig. 4. As seen from Fig. 4a the zero point of the mag- 
nitude difference is at AV = —0.05 mag. The trend 
seen in Fig. 4b indicates a disagreement between our 
values and Maciejewski & Niedzielski's (2007). After 

^ttp:/ /iraf.noao.edu 
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Table 1 Coordinates of the subfields and log of observations with dates, airmass and exposure times for each filter. 
Exposure times are given in seconds. 



Sub- B Band V Band 



field 


«2000 


l$2000 


Exposure 


Date 


Airmass 


Exposure 


Date 


Airmass 


Fl 


04 10 10 


+49 33 20 


30x3 


10/09/2004 


1.036 


10x3 


10/08/2004 


1.257 


Fl 


04 10 10 


+49 33 20 


600x3 


10/09/2004 


1.030 


30x3 


10/08/2004 


1.228 


F2 


04 09 32 


+49 33 20 


30x3 


10/09/2004 


1.025 


10x3 


10/08/2004 


1.188 


F2 


04 09 32 


+49 33 20 


600x3 


10/09/2004 


1.026 


30x3 


10/08/2004 


1.168 


F3 


04 10 10 


+49 27 15 


30x3 


10/09/2004 


1.089 


10x3 


10/08/2004 


1.091 


F3 


04 10 10 


+49 27 15 


600x3 


10/09/2004 


1.063 


30x3 


10/08/2004 


1.079 


F4 


04 09 32 


+49 27 15 


30x3 


10/09/2004 


1.097 


10x3 


10/08/2004 


1.137 


F4 


04 09 32 
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10/09/2004 


1.037 
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1.121 
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Fig. 1 Finding chart of the stars in the field of NGC 1513. Four subfields (Fl, F2, F3 and F4) were also shown. Filled 
circles of different sizes represent brightness of the stars. Smallest size denotes stars of V ~ 19 mag. 



•1 



Table 2 Mean photometric errors and number of stars in 
optical and near-infrared magnitudes and colours for differ- 
ent magnitude intervals. N is the number of the stars for 
each interval. 



V ay (Jb-V 


N w 






N 2 aiass 


(mag) (mag) (mag) 




(mag) 


(mag) (mag) 




(10,11] 0.001 0.001 


1 


0.024 


0.037 


0.033 


1 


(11,12] 0.001 0.001 


7 


0.023 


0.036 


0.031 


7 


(12,13] 0.001 0.001 


9 


0.024 


0.038 


0.033 


9 


(13,14] 0.001 0.001 


32 


0.026 


0.041 


0.036 


32 


(14,15] 0.001 0.001 


58 


0.026 


0.041 


0.037 


58 


(15,16] 0.002 0.003 


70 


0.027 


0.043 


0.040 


70 


(16,17] 0.003 0.005 


105 


0.032 


0.051 


0.051 


105 


(17,18] 0.005 0.009 


133 


0.038 


0.060 


0.067 


133 


(18,19] 0.009 0.016 


194 


0.056 


0.087 


0.109 


191 




Fig. 3 Photometric errors (in magnitude) corresponding 
to the brightness measurement at V and (B — V) are plotted 
against the V-band brightness. Error on the y-axis represent 
the internal errors as estimated by DAOPHOT routine. 
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Fig. 2 Magnitude (a) and colour (b) differences between 
the values calculated in this study and Del Rio & Huesta- 
mendia's (1988) and Frandsen & Arentoft's (1998). V and 
(B — V) represent the magnitudes and colours obtained in 
this study. The dashed lines denote ±lcr. 
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Fig. 4 Magnitude (a) and colour (b) differences between 
the values measured in this study and Maciejewski & 
Niedzielski's (2007) values. V and (B — V) represent the 
magnitudes and colours obtained in this study. The dashed 
line in panel (a) denotes AV = —0.05 mag. 
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analyzing and comparing our data with previous stud- 
ies, we conclude that our values are not in agreement 
with values given by Maciejewski & Niedzielski (2007). 
This discrepancy in Fig. 4 might exist because of the 
fact that Maciejewski & Niedzielski (2007) did not pick 
standard stars from a large magnitude and colour in- 
terval. 

2.2 Near-Infrared Data 

The near-infrared JHK S photometric data were taken 
from the digital Two Micron All-sky SurvcjH (2MASS). 
2MASS uniformly scanned the entire sky in three 
near-infrared bands J (1.25/im), H (1.65/im) and K s 
(2.17/im) with two highly automated 1.3-m telescopes 
equipped with a three-channel camera, where each 
channel consists of a 256 x 256 array of HgCdTe de- 
tectors. The photometric uncertainty of the data is 
less than 0.155 at K s ~ 16.5 magnitude which is 
the photo metric completeness of 2MASS for stars with 



\b\ > 25° (jSkrutskie et all 120061 ). We picked stars in 
a field, of size 25 arcmin 2 , in the direction of NGC 
1513 from Cutri et al.'s Point Source Catalogue (2003) 
and calculated the limiting magnitudes and photomet- 
ric errors are J = 16.5 ± 0.125, H = 16.0 ± 0.143, 
K s = 15.5 ± 0.175 mag. These photometric errors are 
in agreement with the error values given for high lati- 
tude star fields. 

We detected 609 stars in BV photometry in our ob- 
servations. We then compared our results with Cutri 
et al.'s Point Source Catalogue (2003) and matched the 
appropriate stars. After that, we obtained the 2MASS 
magnitudes for the 606 observed stars. 534 out of these 
606 stars have a "AAA" quality flag, which means the 
signal noise ratio is SNR > 10, i.e. they have the 
highest quality measurements. The 2MASS magnitude 
and colour errors are as follows: 0.018 < a,j < 0.145, 
0.032 < <jj-H < 0.232 and 0.030 < (Xj-k b < 0.241 
mag. Errors are given in Table 2 and shown in Fig. 
5. Table 2 reveals that the accuracy of optical data is 
better than near-infrared data, because there are more 
observations in the former. The coordinates, optical 
and near-infrared magnitudes and their errors for 609 
observed stars are given in Table 3. 



3 Data Analysis 

3.1 Cluster's Centre and Radial Density Profile 

NGC 1513 is a cluster with low central concentra- 
tion. The centre of the cluster can only be deter- 
mined roughly by eye- estimation. To determine the 



centre more precisely, we applied the star-count method 
and divided the cluster into one arcmin sized squares. 
We then calculated the surface distributions of those 
squares (Fig. 6). We assumed the centre of the cluster 
as the square's centre with maximum star density. The 
star symbol seen in Fig. 6 indicates the centre of the 
cluster with equatorial coordinates a20oo = 04 /l 09 m 36 s , 
^2000 = 49°28'43" and galactic coordinates I = 152°. 57, 
b = — 1°.64. This result is in agreement with Ma- 
ciejewski & Niedzielski's (2007) a 20 oo = 04 ,l 09 m 46 s , 
<>20oo = 49° 28 28 , but somewhat different than Frolov 
et al.'s (2002) a 20 oo = 04 ,l 10 m 38 s , S 20 oq = 49°3l'00" 
central equatorial coordinates. Considering our results 
were obtained from both optical and near-infrared pho- 
tometric systems, we can state that our central coordi- 
nates are more accurate. 

To establish the radial density profile we counted the 
stars with distances r < 1 arcmin from the centre of the 
cluster. We repeated this process up to r = 15 with 1 
arcmin steps. The next step was to subtract the stars in 
previous areas from the later ones, so that we obtained 
only the amount of the stars within the relevant area, 
not a cumulative count. Finally, we divided the star 
counts in 15 fields by the appropriate areas, i.e. the 
areas of the fields those stars belong to. The density 
uncertainties in each field was calculated using Poisson 
noise statistics. 

Using photometric data wc plotted the radial den- 
sity versus angular separation from the cluster's centre 
in Fig. 6. The figure shows the radial density profile 
from the centre of the cluster to a maximum angular 
separation of 15 arcmin. The density shows a maximum 
at the centre p = 8.5 stars/arcmin 2 and then decreases 
down to p = 5.5 stars/arcmin 2 at 15 arcmin. The de- 
crease becomes asymptotical at r ~ 10 arcmin, after 
that point there are a few cluster stars. To determine 
the structural parameters of the clust er more pre cisely, 
we applied the empirical King model ( King 1966h . The 
King model parameterizes the density function p(r) as: 



P(r) = fbg + 



fa 



1 + (r/r c y 



(1) 



2 http://www.ipac. caltech.edu/2MASS/ 



where fb g , fa and r c are background and central star 
densities and the core radius of the cluster, respectively. 
Fig. 6 reveals the background star density fb g = 5.5 
stars/arcmin 2 . We then compared observational val- 
ues with the ones we obtained from King profile us- 
ing minimum \ 2 statistics. The analysis shows that 
fo = 3.00 ± 0.31 stars/arcmin 2 and r c = 3.98 ± 0.55 
arcmin. The degree of freedom of the analysis (dof) 
is 0.11, while its squared correlation coefficient (R 2 ) is 
0.89. In Fig. 6, the dots stand for the observational 
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Table 3 The coordinates, optical and near-infrared magnitudes and their errors for 609 observed stars. The catalogue 
containing this information can be obtained electronically. 



ID «2000 

(hh:mm:ss) 



02000 B B err 

(dd:mm:ss) (mag) (mag) 



V V err 
(mag) (mag) 



H R err K s K Serr 
(mag) (mag) (mag) (mag) (mag) (mag) 



001 04:09:05.16 

002 04:09:05.17 

003 04:09:05.45 

004 04:09:05.72 

005 04:09:05.87 



49:30:07.62 17.849 0.004 

49:29:11.57 20.053 0.014 

49:32:29.04 19.473 0.010 

49:25:47.56 18.853 0.006 

49:24:24.37 19.812 0.007 



16.563 0.003 

18.510 0.009 

18.108 0.007 

17.446 0.005 

17.553 0.005 



13.859 0.031 

15.487 0.050 

15.384 0.055 

14.780 0.045 

15.483 0.113 



13.440 0.036 

14.897 0.066 

15.065 0.071 

14.402 0.049 

15.040 0.090 



13.257 0.032 

14.727 0.100 

14.627 0.086 

14.201 0.075 

14.826 0.134 



605 04:10:29.69 49:30:37.04 18.652 0.005 16.921 0.003 13.892 0.029 13.158 0.036 12.952 0.030 

606 04:10:29.72 49:27:03.29 19.860 0.010 18.271 0.007 15.137 0.048 14.562 0.057 14.503 0.080 

607 04:10:30.08 49:30:28.02 16.268 0.002 15.115 0.001 13.495 0.027 13.192 0.032 13.134 0.035 

608 04:10:31.19 49:23:41.85 20.263 0.004 16.545 0.003 10.005 0.024 8.608 0.031 8.195 0.023 

609 04:10:32.86 49:26:16.21 17.559 0.003 16.452 0.003 13.824 0.026 13.527 0.032 13.309 0.036 
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Fig. 5 Photometric errors (in magnitude) corresponding 
to the brightness measurement at J, (J — H) and (J — K s ) 
are plotted against the J-band brightness. 



density values, whereas solid line represents the King 
profile. The error bars denote the Poisson errors in ob- 
servations. 



Macieiewski k, Niedzielskil (|2007[ ) analyzed the struc 



tural param eters of 42 open clusters using King model 
( Kingl 1966T ) . They obtained the parameters for NGC 



1513 as: ru m = 9.2, r c = 3.7 arcmin, /o = 2.47 and 
fig = 1.04 stars/arcmin 2 . Comparing our results with 
Maciejewski & Niedzielski's (2007) we can conclude 
that although the limiting radii and the core radii are 
in agreement, the central and background star densities 
somewhat differ. 

To check our background centr al density value, w e 
used the Besangon Galaxy modeH (jRobin et al. 20031 ). 
We assumed the size of the field as 0.05 deg 2 , and the 
magnitude range as 10 < V < 19. We took colour 
excess E(B — V) = 0.68 and distance modulus values 
(V — My) = 12.95 mag from Sections 3.3 and 3.4. Ac- 
cording to Besangon model there are 3.6 field stars per 



in Fig. 6. At the core radius r c = 3.98 arcm in, which w e 
obtained for the cluster using King model 1 Kind 19661 ) . 
the background star density is p = 5.5 stars/arcmin 2 . 
This value differs slightly from Besangon model's 3.6 
stars/arcmin 2 . The observational and model back- 
ground star densities are not in total agreement; there 
is a difference of 1.9 stars/arcmin 2 . This discrepancy 
originates from the galactic model parameters, because 
Galaxy models use several paramet ers with fix e d val - 

~ (I2OO8I ) 



ues for the entire Galaxy. However, iBilir et al 



showed that galactic model parameters are a function 



3 http:/ /model. obs-bcsancon.fr/ 
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of absolute magnitude, galactic latitude and longitude, 
i.e. they do not have fixed values. This means models 
using fixed galactic model parameters can not always 
be used to explain observational data in different direc- 
tions of the Galaxy. 

Consequently, the analyses in this section show that 
the equatorial and galactic coordinates of the centre are 
a 2 ooo = 04' 1 09 m 36 s , 52000 = 49°28'43" and I = 152°. 57, 
b = — 1°.64, respectively, and that the cluster has an 
angular radius of approximately 10 arcmin. 



magnitudes down to the limiting magnitude of V = 16 
mag. The calibration equations used in separating gi- 
ants from dwarfs are as follows: 



Jo = 0.957 x Vb - 1.079, 



H = 0.931 x Vq - 1.240, 



(2) 



(3) 



3.2 Colour-Magnitude Diagrams 

We established optical and near-infrared colour mag- 
nitude diagrams (CMDs) for NGC 1513. In Fig. 7 
we present (B - V)/V, (J - H)/J and (J - K s )/K s 
diagrams for NGC 1513. Since the central concentra- 
tion of the cluster is relatively low, the determination of 
whether if a star is a field star or a member of the clus- 
ter using radial stellar density profile is tough. Even 
though we calculated the core radius of the cluster as 
r c = 3.98 arcmin from Fig. 6, to determine cluster 
members with even more probability we chose the stars 
within the circle with radius r = 5 arcmin. By do- 
ing that, we obtained a more precise main-sequence in 
the CM diagram, which contains 343 stars in our sam- 
ple. All of these 343 main-sequence stars are found 
in-between the two dashed lines in Fig. 7a-c. However, 
there are 110 more stars in-between the dashed lines 
with 5 < r < 15 arcmin. This means the contamina- 
tion in Fig. 7 is about 24 %. 

However, Frolov et al.'s (2002) astrometric study re- 
garding 333 stars in the direction of the cluster provides 
us information about cluster membership of the 609 ob- 
served stars. The filled and open circles in Fig. 7 rep- 
resent Frolov et al.'s (2002) probable cluster members 
[p > 50%) and other observed stars, respectively. The 
reason we use the probable cluster members is to de- 
cide which isochrone to use when determining the age 
of the cluster. The problem with the astrometrically 
determined high probability stars is that they can only 
be detected in brighter magnitudes. As seen in the CM 
diagrams in Fig. 7, there are a few giants in the direc- 
tion of the cluster. The cluster membership probability 
of stars with numbers 360 and 380, which appear in the 
giant regions of the CMDs, are 90% and 93% (Frolov 
et al., 2002) according to their proper motion studies, 
respectively Since there are no spectroscopic studies 
regarding these systems, it is unknown if these stars are 
giants or main-se quence star s . Rec ently, a new method 
was suggested bv lBilir et al.l (|2006l ) to separate field gi- 
ants from field dwarfs. This new method is based on 
the comparison of the 2MASS J, H, K s with the V 



(K s ) = 0.927 x V - 1.292, 



(4) 



where "0" index denotes the de-reddened magnitudes. 
To apply this method, magnitudes of stars should be 
de-reddened. The colour excesses and the de-reddening 
method are given in detail in Section 3.3. The calibra- 
tion and the de-reddened Vb, Jo, Hq, (K s )q magnitudes 
of the two stars were shown in Fig. 8. The de-reddening 
procedure is explained in detail in the next section. As 
seen from the Fig. 8 the stars appear to the right of 
the calibration lines, which means they are in the giant 
region. These two stars are cluster members astromet- 
rically and giant stars photometrically. According to 
Pickles (1998) synthetic data the giant stars numbered 
360 and 380 belong to Kl and K0 spectral types, re- 
spectively 

By selecting the stars within the circle with radius 
r = 5 arcmin and using Frolov et al.'s (2002) high 
probability member stars we determined a more precise 
main-sequence, which is an important factor in decid- 
ing the distance modulus of the cluster. By making this 
selection, we made a precise determination of the clus- 
ter's main-sequence's turn-off point which is the best 
age indicator. 

3.3 Two Colour Diagrams and Colour Excesses 

We present optical and near-infrared two-colour dia- 
grams for r < 5 arcmin in Fig. 9. In Fig. 9a, we 
plotted (J — H)o versus (B — V)o whereas in Fig. 9b 
the axes are (J — K s ) and (B — V) . Namely, we plot- 
ted a near-infrared colour versus an optical colour in 
each panel. To determine the reddening in the dircc- 
ti on of the cluste r we made use of the synthetical library 
of lPicklesI (|l998h . In Pickles' (1998) library we selected 
the metallicity to be [M /H] = dex and main-sequence 
stars of different spectral types and obtained datasets 
for (B - V), (V - J), (V - H), {V - K). Since Pick- 
les' near-infrared data are in Johnson's system, we used 
Carpenter's (2001) transformation equations (Al, A2, 
A3, A4) to obtain magnitudes of 2MASS bands (Table 
4). We plotted the standard main-sequence from Table 



8 




- Besancon 



Q I I i i i I i i i i I i i i i L 

5 10 15 

Radius (arcmin) 

Fig. 6 Radial surface distribution of NGC 1513. Dots: Observational data, solid line: King profile, dashed line: Number 
of background stars calculated using Besangon galaxy model, thick dashed line: Number of observational background stars. 
Errors are determined from sampling statistics (1/y/N where N is the number of stars used in the density estimation at 
that point). The contour map in the upper right shows how we determined the centre of the cluster. The star symbol in 
the upper right represents the centre. 
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Fig. 7 Optical and near-infrared colour-magnitude diagrams for NGC 1513. Empty circles and filled circles denote the 
stars in this study and Frolov et al.'s (2002) high probability cluster stars, respectively. The dashed lines were obtained 
using the cluster members with r < 5 arcmin, while the stars in-between the lines have r < 15 arcmin. 

The stars numbered 360 and 380 are also shown in all panels. 
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Fig. 8 The stars numbered 360 and 380 in magnitude-magnitude diagrams. The dashed lines represent the border of 
giant/dwarf separation, (a) Vo/Jo, (b) Vo/Ho, (c) Vo/(K a )o- 



Table 4 Pickles' (1998) original synthetic colours (1-5). 
(J — H) and (J — K s ) were calculated using Carpenter's 
(2001) transformation equations (6-7). 





Pickles' colours 




2MASS 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


SpTypc 


(B-V) 


(V - J) (V - H) 


(V-K) 


(j — H) (J — K s ) 


05V 


-0.380 


-0.737 


-0.992 


-1.096 


-0.202 


-0.214 


09V 


-0.331 


-0.718 


-0.899 


-0.989 


-0.182 


-0.194 


BOV 


-0.342 


-0.695 


-0.834 


-0.856 


-0.163 


-0.175 


B1V 


-0.244 


-0.632 


-0.732 


-0.780 


-0.143 


-0.145 


B3V 


-0.201 


-0.482 


-0.734 


-0.798 


-0.123 


-0.115 


B5-7V 


-0.139 


-0.322 


-0.370 


-0.395 


-0.094 


-0.076 


B8V 


-0.109 


-0.229 


-0.259 


-0.283 


-0.074 


-0.046 


B9V 


-0.044 


-0.141 


-0.147 


-0.160 


-0.055 


-0.027 


A0V 


0.015 


0.003 


0.000 


-0.009 


-0.045 


-0.017 


A2V 


0.029 


-0.017 


-0.149 


-0.136 


-0.035 


0.003 


A3V 


0.089 


0.036 


-0.129 


-0.126 


-0.016 


0.032 


A5V 


0.153 


0.288 


0.351 


0.360 


0.014 


0.062 


A7V 


0.202 


0.369 


0.459 


0.474 


0.043 


0.101 


FOV 


0.303 


0.533 


0.663 


0.683 


0.082 


0.140 


F2V 


0.395 


0.607 


0.621 


0.645 


0.122 


0.190 


F5V 


0.458 


0.827 


1.030 


1.058 


0.180 


0.248 


F6V 


0.469 


0.900 


1.110 


1.128 


0.210 


0.288 


F8V 


0.542 


1.012 


1.219 


1.245 


0.249 


0.317 


GOV 


0.571 


1.017 


1.275 


1.300 


0.298 


0.376 


G5V 


0.686 


1.190 


1.452 


1.507 


0.288 


0.386 


K2V 


0.924 


1.650 


2.142 


2.226 


0.445 


0.563 


K3V 


0.930 


1.808 


2.366 


2.396 


0.484 


0.612 


K4V 


1.085 


1.973 


2.541 


2.643 


0.523 


0.661 


K5V 


1.205 


2.172 


2.770 


2.872 


0.553 


0.691 


K7V 


1.368 


2.398 


2.968 


3.094 


0.602 


0.779 


MOV 


1.321 


2.855 


3.514 


3.678 


0.612 


0.809 


M1V 


1.375 


2.974 


3.622 


3.896 


0.602 


0.909 


M2V 


1.436 


3.294 


3.961 


4.139 


0.602 


0.829 


M3V 


1.515 


3.817 


4.449 


4.674 


0.582 


0.839 


M4V 


1.594 


4.422 


5.038 


5.306 


0.563 


0.860 


M5V 


1.663 


5.253 


6.096 


6.393 


0.563 


0.920 


M6V 


1.816 


6.362 


7.026 


7.409 


0.602 


1.008 



4 and our observations in two-colour diagrams (Fig. 9) . 
In order to determine the reddening we calculated loci 
(represented by star symbols in Fig. 9) for our obser- 
vational data and plot the best fit for those loci. These 
loci represent the main-sequence for our observations 
and we slide that main-sequence in both directions un- 
til it fits best with the standard main-sequence. The 
amount of sliding gives us the colour excesses for Fig. 
9a and b, E(J-H), E(B-V) and E(J-K S ), E(B-V) 
respectively. The colour excesses and their relative er- 
rors we obtained using minimum \ 2 method arc as fol- 
lows: E(J -H) = 0.21 ± 0.02, E(B - V) = 0.68 ± 0.06 
and E(J - K s ) = 0.33 ± 0.04, E(B — V) = 0.68 ± 0.06 
mag, respectively. The confidence level of colour excess 
errors is 99.5%. The contour maps of two-colour dia- 
grams in Fig. 10 show the optimum colour excesses for 
B V and 2MASS photometri es. 



Del Rio fe Huestamendial (|l9881 ) obtained RGU colour| 
excess as E(G — R) = 0.94 mag from standard 



stars. iFrolov et al 
E(B - V) 
culated E(B - V) 



(2002) converted this value into 
using Steinlin's (1968) formula and cal- 
0.67 mag. The recent study 



of iMacieiewski fe Niedzielskil |2007l) estimates E(B 
V) = 0.76to;ig mag. Obviously, our value is almost in 
perfect agreement with Frolov et al.'s (2002), whereas 
Maciejewski & Niedzielski's (2007) agrees only within 

er ror bars. 

Fiorucci fc Munari ( 2003[ ) calculated the colour ex- 



cess values for 2MASS photometric system. According 
to them the colour excess ratios are: E(J — H)/E(B — 
V) = 0.322 ± 0.074 and E(J - K S )/E(B - V) = 
0.505 ± 0.053 mag. Wc ended up with the following 
results: E{J - H)/E(B - V) = 0.309 ± 0.130 and 
E(J - K a )/E(B - V) = 0.485 ± 0.150 mag. The 
results produced for normal interstellar medium by 
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Fig. 9 Optical and near-infrared two-colour diagrams. Filled circles, the solid line and star symbols represent stars within 
r < 5 arcmin angular separation from the cluster centre, Pickles' (1998) original synthetic main-sequences and the loci of 
the main-sequence stars for the cluster, respectively. 




Fig. 10 Colour excess diagrams. The plus symbol represents the reduced minimum \ 2 m both panels, whereas the 
contours denote the probability distribution of reduced x 2 . These contours are used in determining the colour excesses and 
their relative errors. 
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Fig. 11 V/(B - V), J/(J - H) and J/ (J - K„) CMDs of the cluster. The black circles represent high probability 
(P > 50%) cluster members; whereas the grey circles denote the low probability cluster members and field stars. Thin solid 
line: ZAMS, thick solid curve line: \og(t/yr) = 8.40, thin dashed curve line: \og(t/yr) = 8.60. 
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Fig. 12 V/(B — V), J/(J — H) and J/(J — K a ) CMDs of the cluster. The solid, dashed upper and lower curves correspond 
to \og(t/yr) — 8.40, \og(t/yr) = 8.35 and \og(t/yr) = 8.44 isochrones, respectively. Finally, the thick-dashed line stands for 
binary sequence. 



Fiorucci fe Munari ( 20031 ) is in agreement with our re- 
sults. 

3.4 Distance and Age of the Cluster 

The zero age main-sequence (ZAMS) fitting procedure 
was used to derive the distance to the cluster. We added 
Schmidt-Kaler's (1982) ZAMS to the optical CMDs in 
Fig. 11. We used the colour excess E{B — V) = 0.68 
mag discussed in Section 3.3 and reddened Schmidt- 
Kaler's (1982) ZAMS accordingly. We slid the cluster's 
main-sequence vertic ally until it overlapped with the 
Sch midt- Kalerl (|1982l ) one. The distance modulus of 
the cluster is the amount of sliding we have applied to 
the cluster's main-sequence, which is (V — My) = 12.95 
mag. Since there is no 2MASS data regarding Schmidt- 
Kaler's (1982) ZAMS, we have selected proper 2MASS 
data from the Padova isochrone with solar metallic- 
ity. We plotted 2MASS ZAMS dMarigo et all l2008h 
in Fig. 11. Then, we applied the previous procedure 
to 2MASS CMDs (Fig. 11) using reddening values of 
E(J-H) = 0.21 and E{J -K s ) = 0.33 mag. The opti- 
cal and near-infrared ZAMS are given as thin solid lines 



4 http://stev.oapd.inaf.it/cgi-bin/cmd 



in all panels of Fig. 1 1 . Finally, we calculated the dis- 
tance moduli for J/(J ~ H) and K S /(J — K s ) diagrams 
and obtained (J - Mj) = 11.40 and (K s - M Ks ) = 
10.95 mag, respectively. To de-redden the distance 
moduli calculated from CMDs, we used Fiorucci & 
Munari's (2003) formulae: A v = 3.1 x E(B - V), 
Aj = 0.887 x E(B - V), A Ks = 0.322 x E(B - V) 
and obtained: (V - M v )o = 10.84 ±0.19, (J - Mj) = 
10.80 ± 0.10 and (K s - M K Jo = 10.73 ± 0.10 mag, 
using these de-reddened distance moduli we obtained 
the distance of NGC 1513 1472, 1445 and 1400 pc, 
respectively. These results, obtained from two differ- 
ent photometric systems, are of ±5% re lative differ 
ence with each other. Moreover, both Frolov et al 



( 2002h and lMaciejewski fc Niedzielskil (|2007h calculated 



the distance of the cluster as 1320 pc, which indicates 
a relative difference of ~ 10%. 

The age of a star cluster can be determined by com- 
paring the observed CMDs with theoretical isochrones. 
To determine the age of a cluster colour excess, dis- 
tance modulus and metallicity needs to be known. In 
this study, we determined the colour excess and dis- 
tance modulus, separately. However, since there is no 
spectroscopic data regarding the NGC 1513 its metallic- 
ity remains unknown. Hence, we took into account the 
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mctallicity during the calculation of the age by using 
a set of isochroncs for stars with masses 0.15 < Mq < 
100, different metal abundance 0.0001 < Z < 0.03 and 
ages from \og(t/yr) < 10.24 published on the web site 
of the Padova research group and described in the work 
of Marigo et alj ( 20081 ). We produced several isochrones 
for optical and near-infrared bands with different metal 
abundances ranging from 0.0095 to 0.024, which cor- 
responds to [M/H] = -0.30 and [M/H] = 0.10 dex, 
respectively (solar abundance was assumed as Z Q = 
0.019). The CMDs and the isochrones are given in Fig. 
11. On each CMD, we plotted three isochrones repre- 
senting three different ages: ZAMS, \og(t/yr) = 8.40, 
\og(t/yr) — 8.60. As seen from Fig. 11, the isochrones 
with Z = 0.015 provide the best fit for the cluster's 
main-sequence, main-sequence turn-off point and giant 
stars. Therefore, we assumed Z = 0.015 (corresponds 
to a mctallicity of [M/H] = —0.10 dex) to be the metal 
abundance of the cluster. 

To determine the age precisely, we plotted three 
isochrones with ages \og{t/yr) = 8.35, log(t/yr) = 8.40 
and \og(t/yr) = 8.44 in the CMDs with Z = 0.015 
(Fig. 12). The isochrone with \og(t/yr) = 8.40 seemed 
to represent the cluster best, because it fits with both 
the main-sequence and the giants of the cluster. The 
loop (Fig. 12a: B — V = 0.61, V = 11.93; Fig. 12b: 
J - H = 0.16, J = 10.44; Fig. 12c: J — K s = 0.26, 
K s = 10.05 mag) in the \og(t/yr) = 8.40 isochrone 
marks the ending of hydrogen burning in the core, com- 
pression of the core and hydrogen burning in the thick 
layer. The star numbered 380 is located almost at the 
base of the red giant branch, and the star numbered 
360 at the stage of helium burning. The turn-off point 
of the main-sequence corresponding to this isochrone is 
(B — V)o = —0.11 mag, which is the colour index for 
the B8 spectral type. To determine the binary effect 
regarding the cluster, we assumed the stars in the clus- 
ter to be binary systems of equal-massed components. 
In this case, absolute magnitude is 0.75 mag brighter 
than it should normally be. This binary main-sequence 
with \og(t/yr) = 8.40 has been plotted in Fig. 12 as a 
thick-dashed line. 

Consequently, the average distance of the cluster is 
d = 1440 ± 80 pc, the mctallicity Z = 0.015 ± 0.004 
([M/H] = -0.10 ± 0.10 dex) and average age is 
\og(t/yr) ~ 8.40. These results partially agree with 
Frolov et al.'s (2002) results. They determined the dis- 
tance, mctallicity and average age as d = 1320 pc, Z = 
0.019 ([M/H] = dex) and \og(t/yr) = 8.40. Com- 
paring our results with Frolov et al.'s (2002), we can 
claim that the distance we calculated is within 10% of 
each other , while the metallicity seem s lower and age is 
the same. Macieiewski fc Niedzielskil (|2007t ) estimated 



the age of the cluster \og(t/yr) = 7.4, which is not in 
agreement with the age we calculated, which might be 
due to the fact that they did not take the giant stars of 
the cluster into account. To determine the integrated 
absolute magnitudes and colours of NGC 1513 in opti- 
cal and near- infrared, we used Lata et al.'s (2002) equa- 
tions and obtained the following results for cluster stars 
(r < 5 arcmin): I(M V ) = -2.136, I(B - V) = 0.032, 
I(J - H) = 0.053 and I(J - K s ) = -0.037 mag. 
According to Table 4 the optical integrated colour cor- 
responds to A2 spectral type. 



4 Conclusion 

We present CCD BV and JHK S 2MASS photometric 
data for the low central concentration young star cluster 
NGC 1513. The results obtained in the analysis are the 
following: 

i) We determined the centre of the cluster as c^ooo = 
04 /l 09 m 36 s , <5 2 ooo = 49°28'43" and its galactic coordi- 
nates I = 152°. 57, b = -1°.64. The radial density 
profile shows that the angular radius of the cluster is 
r = 10 arcmin. 

ii) The optical and near-infrared colours of the clus- 
ter main-sequence reveal the colour excesses, E(B — 
V) = 0.68 ± 0.06, E(J - H) = 0.21 ± 0.02 and 
E(J - K s ) = 0.33 ± 0.04 mag. We estimated E(J - 
H)/E(B-V) = 0.309 ±0.130 and E(J - K S )/E(B - 
V) = 0.485 ±0.150 mag. 

iii) We compared the optical main sequence of the 
cluster with Schmidt-Kaler's (1982) ZAMS the near- 
infrar ed one with Padova isochrones' ZAMS ( Marigo et al.^ 
20081 ). We obtained the distance moduli as (V — 
M v ) = 10.84±0.19 for optical colour and (J-Mj) = 
10.80 ±0.10, (K a -M Ks ) = 10.73 ±0.10 mag for near- 
infrared colours. The average distance derived from 
these moduli is 1440 ± 80 pc. 

iv) The Padova isochrone with Z = 0.015 and 
\og(t/yr) = 8.40 provides the best fit for the cluster 
in both optical and near-infrared CMDs. Therefore, we 
conclude that the cluster is \og(t/yr) = 8.40 ±0.04 and 
has a metallicity of [M/H] = -0.10 ± 0.10 dex. 
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